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ABSTRACT 

Using the cosmic virial theorem, Press-Schechter analysis and numerical 
simulations, we compute the expected X-ray background (XRB) from the diffuse 
IGM with the clumping factor expected from gravitational shock heating. 

The predicted fluxes and temperatures are excluded from the observed XRB. 

The predicted clumping can be reduced by entropy injection. The required 
energy is computed from the two-point correlation function, as well as from 
Press-Schechter formalisms. The minimal energy injection of 1 keV/nucleon 
excludes radiative or gravitational heating as a primary energy source. We argue 
that the intergalactic medium (IGM) must have been heated through violent 
processes such as massive supernova bursts. If the heating proceeded through 
supernova explosions, it likely proceeded in bursts which may be observable in 
high redshift supernova searches. Within our model we reproduce the observed 
cluster luminosity-temperature relation with energy injection of 1 keV/nucleon 
if this injection is assumed to be uncorrelated with the local density. These 
parameters predict that the diffuse IGM soft XRB has a temperature of ~ 1 
keV with a flux near 10 keV/cm^ s str keV, which may be detectable in the near 
future. 

Subject Headings: cosmology: theory and diffuse radiation - X-rays: general - 
supernova: general 


1. Introduction 

The intergalactic medium has been a mysterious subject for many decades. While the 
baryonic content of the universe is still uncertain, the current Big-Bang nucleosynthesis 
value of VLi)h? = 0.02 (Buries and Tytler 1998) is much larger than the observed stellar 
density = 0.0025 (Fukogita et al. 1997). In recent years, the large gas mass in clusters 
of galaxies (White et al. 1993) has uncovered a diffuse reservoir of baryons. Since the gas 
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fraction in clusters should be representative of the global mean, this hot diffuse component 
of baryons could account for the majority of baryons in the universe. Barcons et al. (1991) 
hrst applied the soft X-ray background (XRB) constraints to the IGM. Since their work, 
the major portion of the XRB has been resolved into point sources (Basinger et al. 1993), 
strengthening the bound. We here present new aspects to the argument by estimating the 
temperature with the cosmic virial theorem, and calculating the clumping factor directly 
through numerical simulations and Press-Schechter estimates. These two ingredients 
rule out a passively evolving IGM where the energetics are dominated by gravity. We 
discuss these constraints in more detail in section Using our galaxy as a representative 
gravitational well, we place constraints on the non-gravitational heat injection required to 
explain the small observed XRB contribution from our halo in section In section ^ we 
discuss predictions and implications of the model. It explains the luminosity-temperature 
relation in clusters of galaxies, and predicts a very sudden heating period for each galaxy. 


2. Cosmological X-ray background 


The cosmological background flux is 


,7 n2 2 l<^T'/^P){l + zf .^dr . 

f = ^ - d,-dz 


erg 


cm^ s str 


( 1 ) 


where cIa and denote the angular diameter and luminosity distance respectively, related 
by di = {1 + z^dA- Pc = SHq/StiG is the critical closure density of the universe. 6 is 
the local overdensity of the gas which is related to (d^) = ^(0) in terms of the correlation 
function at zero lag. The volume emission coefficient esyoi ~ 1-7 x 10“^^cz/0.875mp ([Rybicki 
Ik. Lightman 19791) is expressed in terms of the proton mass rup at a solar hydrogen-helium 
mixture with a coefficient cz ~ Z/ZQ{A'keV/T) -|- 1 to correct for metal cooling (Raymond 
et al. 1976). The approximation for cz is accurate to 20% in the interval 0.01 keV < T < 
10 keV. We have assumed a mean gaunt factor gs ~ 1.2. We will use the Z = Z^jA as the 
lower plausible limit on the intergalactic medium inferred from clusters of galaxies. Higher 
metallicities strengthen the bounds derived here. 

In a flat universe with a powerlaw correlation function ^ = (r/ro)“^'®/(l -|- zY the 
characteristic temperature evolves as oc (1 -|- and is expected to evolve more 

slowly in a low density model. To simplify (|l|), we neglect this weak redshift dependence. 
We express the integrand in terms of tq = to77o) Ihe age of the universe in units of the 
Hubble constant. Recall that Tq = 2/3 for a flat universe and tq = 1 for an empty universe. 
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Since F = i{Q)nlple^TQc/H q, we find in differential energy units 


F = m 



T 


0.13 keV 


- 1/2 


If 

h 


exp( 


/iz/ /keV \ 0.16 keV 

kT y T y cm^ s str keV 


( 2 ) 


X-ray emission is proportional to the square of the density. Its emission weighted 
temperature is well approximated by the pair weighted temperature given be the cosmic 
virial theorem (Peebles 1980) (locfkpc) ^ 0.128 keVwhere is a short scale 

cutoff of the gas autocorrelation function. From the cluster abundance (Pen 1998b), we use 
(Tg = .53f2“-^®, and use the approximation {T)x = 0.13 keV. 


The observed extragalactic XRB at 1/4 keV has an upper limit of 15 keV/cm^ s keV 
str (Cui et al. 1996). From shadowing studies, they found an extragalactic background 
of 45 keV/cm^ s keV str, of which at least 30 keV/cm^ s keV str has been resolved into 
point sources (Hasinger et al. 1993). Using To/h = 1 and taking into account a factor of 2 
suppression of the flux due to the energy dependent Gaunt factor (Rybicki and Lightman 
1979), we obtain the main constraint 


e(o) < 60 . 


( 3 ) 


A less certain prediction is the clumping factor ^(0). For a scale invariant spectrum in a 
flat universe, this is expected to be constant in time. According to Press-Schechter theory, 
most of the mass is in collapsed and virialized objects of approximately the non-linear 
mass scale. The average overdensity of collapsed objects is at least 178, and generally 
higher if the objects are not homogeneous or have collapsed at higher redshift. We expect 
at least a factor of two higher clumping, since virialized objects are certainly not uniform 
density objects. Estimates using stable clustering (Jain 1997) and numerical pure N-body 
simulations indicate that at least for the dark matter component, ^(0) 10^. With 

these values, we far exceed the observed diffuse soft XRB constraints. If the distribution 
of gas traces light down to scale r^, and then forms a constant density core, we hnd 
^(0) = (rg/5.5h“^Mpc)“^'®. Using (|) we hnd Vg > 800h“^ kpc, i.e. gas does not cluster on 
scales shorter than a Mpc, even though all other forms of matter do. 

We can compare these estimates with direct numerical simulations (Cen et al. 1995). 
Scaling their Figures 3 and 4 to the cosmological XRB by multiplying the vertical abscissa 
by 5 X 10^® to obtain units of keV/cm^ s str keV, we hnd that they predict emission just 
below the observed upper bounds. Emission is weighted by the smallest scale clumping, 
which simulations at hnite resolution always underestimate. We can see the ehect in their 
Figure 8. Their COBE normalized simulation should predict weighted temperatures 
of ~ 0.3 keV, while the simulated temperatures are signhcantly higher. Increasing the 
resolution would primarily enhance the emission of smaller mass and lower temperature 
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objects, which are smaller and more subject to numerical limitations. It would also raise 
the clumping factor and X-ray emission. More recent simulations may move the results in 
the expected direction (Cen and Ostriker 1998). 

Using the cosmological Moving Mesh Hydro Code (Pen 1998a), we have performed 
direct simulations to measure ■^(0)- We use a scale-free CDM models with 128^ grid cells 
and an equal number of dark matter particles. The correlation function is oc (r/ro)“^'® 
with ro = 4 h“^Mpc in a 16 h“^Mpc box. In order to extend the dynamic range, the initial 
conditions were determined by factoring the correlation function (Pen 1997). The grid was 
allowed to compress a factor of 10, giving a minimal grid spacing of 13h~^ kpc, with an 
effective resolution which is several times that value. We used Qb = 0.1. In Figure we plot 
the gas correlation function ,^(r) at different redshifts scaled to their relative amplitudes 
using scale invariance. The simulation moves from the lower right to the upper left as 
the non-linear mass scale grows, and the effective resolution increases. We see that the 
clumping factor is limited by the initial resolution at all times. The clumping factor ,^(0) 
monotonically increases with time, reaching a value of ~ 900 at z = 0 when the non-linear 
mass scale reaches 32 grid cells. Our limited resolution allows us to only place a lower 
bound on ^(0) 900. One would expect low density models to have higher clumping factors 

since objects formed earlier. 

We must thus conclude that the gas could not have followed the gravitational evolution 
of the dark matter, as one would have expected in the scenario of purely gravitationally 
driven collapse. Two potential solutions come to mind. Baryons could have cooled into 
compact objects which do not contribute to the XRB. But this poses problems for clusters 
of galaxies, where vast quantities of gas have been observed (White et al. 1993). In a 
hierarchical model of structure formation, clusters form from the merger of smaller objects. 
But if the smaller objects had their baryons locked up into compact objects, it would be 
very difficult to release the gas again to form the Intra-cluster medium (ICM). 

The second solution would be to heat the gas sufficiently, such that it does not fall 
into the gravitational potential well of the dark matter halos. The required energies can 
be estimated from the two-point correlation function. We calculate the energy required 
to convolve the gas with a spherical tophat to satisfy (^. For a gas distribution pg in a 
dark matter potential 0, we have E = J pgcj), which for our power law correlation formally 
diverges due to long wavelength contributions. We consider a specihc model in which a 
smoothed density distribution is given as a convolution over the original gas held. The 
change in energy is hnite, in fact AE = f Apg(j) and the average specihc energy change is 
AS = HoPcdvrGrQ / hF(r/ro)^(r/ro)d^(r/ro). We let IF be a spherical tophat of radius Vg. 
Using as from above, we obtain AS = 0.8 x 3(5 — n)“^(3 — n)“^(2 — n)“^(rg/ro)^“” keV, 
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where n ~ 1.8. The constraint for the smoothed gas held simply translates to ^ 60, 
which can again be expressed in terms of the correlation fnnction Vg ^ kpc. 

In this model, approximately 2 keV of heating wonld be reqnired to satisfy the constraints 
from the XRB. This assnmes that smoothing was applied to all points in space, which 
probably overestimates the energy reqnirement. We obtain a lower energy bonnd by jnst 
trnncating the correlation fnnction, so that ^(r) = (rg/ro)“^'® for r < Vg. This reqnires 0.25 
keV, bnt mnst be taken as a lower limit since it violates locality. We will now proceed to 
bnild a more realistic model which allows ns to qnantify the heat injection reqnired to evade 
the bonnd from eqnation (||). 


3. Halo Model 


We will consider a model for the distribntion of gas within each halo. We assnme 
the gravitational potential to be dominated by dark matter with a singnlar isothermal 
sphere mass distribntion. The mass-rotation relation for the dark matter is taken from the 
isothermal Press-Schechter model (Eke, Cole and Frenk 1997) M = ■ '^c is 

the circnlar rotation speed of the object. The virial radins ta of an object is dehned as the 
hdncial point where the mean interior density is A ~ 17812“^/^ times the mean cosmological 
density. It corresponds to half the tnrnaronnd radins of a tophat model. For an isothermal 
sphere ta = The gas is assnmed to have two phases: an onter region r > ri 

where gas traces mass isothermally with a gas fraction fg, and an inner region where the gas 
has been heated to constant entropy. Snch a distribntion is the generic ontcome of central 
non-gravitational heating where entropy is injected in the center, which is then convectively 
transported ontward. Observationally (Danos and Pen 1998, Cooray 1998), the best £t for 
the gas fraction is fg = 0.0Qh~^^^. For the onter region, the gas density prohle is given by 
an isothermal sphere Pg = fgV^/ATir'^G. In the inner region r < ri, the gas density is then 


1 3/2 

Pg = fffic ^ M • We can dehne an eqnivalent core radins Vc = 0.487ri. For this 

choice, the free-free Inminosity of the two phase object is identical to an isothermal (3 = 2/3 
model with core radins Tc (Jones and Forman 1984). Lx = esfgV/T^^‘^/{16rcG‘^).Dne to 
the enhanced central entropy, the emission weighted temperatnre of onr objects will be 
10% higher than the isothermal temperature {T)x = 1.1 x pmpv//2k.We can estimate 
the energy required to increase ri. We calculate the difference in binding energy between 
a singular isothermal gas distribution in the dark matter well, to one with hnite ri. 
While the energy injection in our model is concentrated to the center, we dehne a mean 
energy 6T to be t he total injec ted energy divided by the total virial gas mass. We obtain 
^ = 1.582 (^1 — — 1.029^^ In the limit of 5T/T <C 1, 


we hnd ^ = 0.814f^.The 
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clumping factor C of an individual object is C = ^ (|l.553^ — l^and in the limit that 
Ti -C T, we have the linear relation C = 0.52TA/5T. We hnd the average clumping factor 
(C) = ^(0) = ^C(T)dT where the lower cutoff assumes that objects with energy- 
injection larger than their virial temperature have ejected all their gas. The Press-Schechter 
distribution function is /(> M) = We use the standard value 6c = 1.686. 

' CT 

Using a power law mass fluctuation spectrum a{M) = (M/Mg)”" (Pen 1998b), we hnd 


{C) = 





--1 


2e '^dw 


( 4 ) 


/Wi 


For most CDM-like spectra, 
+ -\/Ferfc(i/uh') ^where 


where Tg ~ (Pen 1998b) and wi = 6l{6T/Ts)' 
a ~ 1/3, which simplihes the (Q) to (C) = x 

now Wi = 5Q6T/T. The leading order expansion of the bracket in is We 

can compare the predicted clumping to the constraint (||). By assumption, all gas has a 
temperature of ST. For an = 1 scenario, we can satisfy the diffuse XRB constraint at 
1 keV for 6T ^ 1 keV. We used an extragalactic diffuse hux limit of 8 keV/cm^ s str keV 
(Basinger et al. 1993) at 1 keV. 


4. Galactic X-ray background 


While the cosmological constraint depends on parameters including geometry of 
the universe and the evolution of the clumping, our own Galaxy provides a completely 
independent constraint on the thermal state of the intergalactic medium. Simulations 
(Navarro et al. 1995) and observations indicate that the matter distribution is well 
approximated by an isothermal sphere extending from ca 1/lOth of the virial radius out 
to the virial radius. In that range, gas and dark matter do appear to trace each other to 
within a factor of two, and are within a factor of two of the singular isothermal sphere. 
Approximating our solar system to be at the center of the gravitational galaxy halo, we hnd 


F = eT^/Vi 


AnGrl) [3 Jo 25 ^ 


erg 


cm^ s str 


( 5 ) 


To support the outer gas in pressure equilibrium, the central gas will have a temperature 
higher than the virial halo temperature. The emission weighted temperature seen at the 
center differs from the virial temperature, becomes instead Tx = 2.25 x fimpv‘^/2k, more 
than twice the isothermal value. For a rotation speed of 220 km/s, we expect Tx = 0.34keV, 
independent of the core radius ri as long it is greater than the heliocentric distance: 
ri > Rq. 
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The 1/4 keV galactic background is known to be dominated by the local hot ISM 
bubble, which is difficult to disentangle from a halo distribution. Instead, we will use the 
1 keV background constraint, for which (^) yields F = 0.003/(h^r^) keV/cm^ s str keV. ri 
is measured in Mpc. Current reports he near 20 keV/cm^ s str keV (Miyaji et al. 1998), 
of which Hasinger et al. (1993) constrain the homogeneous component to less than 20%. 
This value is lower than the previously used value since the cosmological XRB could be 
clumped. The lower limit for the heated core radius is ri > 150h“^kpc. The virial radius 
of our galaxy is ta = 127h“^ kpc, which is smaller than the required core radius. This 
corresponds to a scenario where all the IGM gas must have been completely expelled from 
our galaxy. Unbinding the gas from the galactic potential well requires a minimum heat 
injection of 150 eV per nucleon. We now have two independent lines of evidence that very 
signihcant non-gravitational heating of the IGM must have occurred. Supernovae might be 
one potential energy source which we will discuss further in section |^. 

In the calculation of our energy injection estimates, we have taken the binding energy 
difference between a model where the gas traces the mass to one where it satishes the 
XRB constraint. It would seem plausible that this is a lower bound to the actual energy 
injection required. In addition to the potential energy change, a similar amount of energy 
may be expended to raise the entropy of the gas. Placing the heating process too early 
causes some of the injected energy to be lost to cosmic expansion, thus requiring even larger 
energy inputs. We should also note that in a cosmological context, energy is not strictly 
conserved, and it is in principle possible to extract energy from the dark matter through its 
time dependent potential. Simulations seem to suggest that gas traces the total mass up 
to overdensities of at least 10^ (Navarro et al. 1995), suggesting that in fact this energy 
exchange process does not play a dominant role. We assumed that objects formed recently. 
An earlier formation time would increase the clumping factor, and raise the predicted X-ray 
emission. 


5. Potential Tests 

We now consider several predictions of the violent heating model. We have assumed an 
equal specihc energy injection in all objects independent of mass. This model predicts the 
core radius Vc oc T~^, which can be compared to the cluster luminosity-temperature relation. 
From the equations in section L 44 = (p^v) (^) ^ lO^^erg/swhich can 

be compared to the Henry and Arnaud (1991) data. Their relation is consistent with the 
slope L (xT^, and agrees for values of 5T ~ 1 keV. Metzler and Evrard (1994) had studied 
numerical simulations with similar parameters and arrived at consistent conclusions. In 



this simplest model where the heating is uncorrelated with the mass of the object, we hnd 
that the required heating to explain the cluster core radii can simultaneously explain the 
soft XRB. 

If the current spatial stellar distribution is a tracer of the heat source distribution, we 
obtain an upper limit on the heating time scale to be the cooling time of the gas. The 
observed metalicities in clusters is about 1/3 solar, suggesting that approximately 1% of the 
gas had undergone nuclear burning, releasing a few MeV of kinetic energy for each burnt 
nucleon. We have a mean energy injection about 10 keV for the gas, of which 10% must 
remain after radiative and adiabatic losses. We require that the heating time not be more 
than 10 times longer than the cooling time in order for the injected heat to be retained. 
The local luminosity density in the disk is ~ 0.07 Lq/pc^ (Binney and Tremaine 1987). 
Our galaxy has a virial mass of 2.6 x IO^^Mq, and using the cosmological gas fraction fg 
from above, we expect a total gas mass of 1.6 x At a total luminosity of 

1.4 X lO^'^L©, we have a mean gas to light ratio of We expect a cooling 

time of fcooi ~ 2.8 x 10^ ^ years assuming a temperature T = 1 keV. Allowing 90% 

of the injected energy to dissipate radiatively, we need to inject the energy in St = lOtcooi- 
We infer that the heating must have occurred over ~ 10^ years. 

Each supernove injects about ten Mq of metals, so we expect ~ 10® supernovae in the 
cooling time. We draw two conclusions: 1. there must have been a short burst of very active 
star formation and supernova explosions, during which about 10 supernovae detonated per 
year. 2. the galaxy would have been at least ten times brighter during this period since 
only stars of M/L<1 will burn out. The epoch of this explosive heating is not constrained, 
and different galaxies may burst at different times. We do predict that a small fraction 
(~ 10“®) of all galaxies at high redshift {z > 1) undergo enormous supernova bursts. Dust 
obscuration can weaken these luminosities. 


6. Conclusion 

We have shown that the current constraints from the XRB signihcantly constrain the 
current thermal state of the intergalactic medium. The absence of signihcant emission from 
our galactic halo places a lower bound of 0.15 keV energy injection. A stronger constraint 
arises from the absence of diffuse intergalactic emission, for which the emission in enhanced 
by the clumping factor of the gas. We found that a mean energy injection of ST ^ 1 
keV is required. Photoheating is not a possible solution. Such a process would raise the 
temperature of the gas to ~ 0.002 keV. In the best possible case, this heating occured at 
~ 0, and the objects instantaneously collapse today. The XRB constraints require that 
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the overdensity be less than 100, for which adiabatic compression can raise the energy 
to at most 0.04 keV, far short of the required heat input. Our halo model also accounts 
for the observed cluster luminosity-temperature relation if 6T ~ 1 keV. The cosmological 
diffuse XRB should be detectable soon at temperatures near 1 keV. This is consistent with 
past simulations (Metzler and Evrard 1994, Suginohara and Ostriker 1998). If the heating 
process traced the present day light distribution, we obtain upper limits on the heating 
time scale, which predicts galaxies undergoing short violent supernova bursts at ^ 1. 
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scale invariant correlation function 



Fig. 1.— Numerical estimates of the correlation function clumping factor {C) = t(0) for 
gas in adiabatic scale free simulations. The solid lines show the correlation function in the 
simulation at redshifts = 180,16,4,1,0 scaled to the appropriate non-linear mass scale. 
The dashed line is the scale free correlation for n = —1.8 and tq = 4 /i“^Mpc. We see that 
the simulations are still limited by resolution, and give a lower bound of {C) ^ 900. This 
rules out passively evolving IGM models, and requires signihcant active heating to evade 
soft XRB constraints. 




